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Abstract 

We investigate instantons in finite temperature QCD via Witten’s holographic 
QCD. To study the deconfinement phase, we use the setup proposed in [T]. We find 
that the sizes of the instantons are stabilized at certain values both in the confinement 
and deconfinement phases. This agrees with the numerical result in the lattice gauge 
theory. Besides we find that the gravity duals of the large and small instantons in the 
deconfinement phase have different topologies. We also argue that the fluctuation of the 
topological charges is large in confinement phase while it is exponentially suppressed in 
deconfinement phase, and a continuous transition occurs at the Gross-Witten-Wadia 
(GWW) point. It would be difficult to observe the counterpart of this transition in 
lattice QCD, since the GWW point in QCD may stay at an unstable branch. 
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1 Introduction 


Although instantons are essential ingredients in QCD, it is difficult to understand their dy¬ 
namics because of the strong coupling nature of the theory. Perturbative calculations are 
justified only for small instantons or at high temperature, and a suitable effective theory 
which describes instantons is not known. Hence we may have to rely on numerical calcula¬ 
tions in lattice gauge theory to illuminate their properties. 

One possible tool for analyzing instantons is holographic QCD proposed by Witten [2]. 
Although holographic QCD is different from real QCD in quantitive details, it has success¬ 
fully explained various qualitative aspects of large-TV QCD [31IH El [6] . (See also [3 [8] for 
recent developments.) Hence we expect that holography can also reveal the natures of in¬ 
stantons. In particular, we focus on the dynamics of instantons at finite temperature in this 
study. 

Through holographic QCD, instantons at low temperature (in confinement phase) have 
been studied in iiiiniin]. It was shown that the energy of an instanton with a particular 
size approaches to zero in the large-limit, which indicates the large fluctuations of the 
topological charge in the confinement phase. This is consistent with the previous theoretical 
insights [12] and numerical calculations in lattice gauge theory [T3] . 

However, dynamics of instantons in deconfinement phase is less clear. Although the per¬ 
turbative calculations work in certain circumstances [ 13 ], the whole instanton dynamics has 
not been understood. Their dynamics around the critical temperature would be particularly 
important to reveal the mechanism of the phase transition, and hence it is interesting to 
study it in holographic QCD. In this direction, the black D4-brane geometry [ 3 ], which was 
supposed to be the gravity dual of the deconfinement phase, has been studied initially. In 
particular, some agreements with the expected properties of the instantons were reported in 
Ref. m- However there are also some disagreements. For example, the instanton density 
n{p,T), which is the vacuum expectation value (vev) of the single QCD instanton with a 
size p at temperature T, shows unexpected behaviours. In holographic QCD, the instanton 
density is calculated from the DBI action of a DO-brane iiiiniiii], and the result in the 
black D4-brane background is given by Cl 



(p<i/r) 
(p ^ i/r) 


( 1 . 1 ) 


where \ym is a dimensionless ’t Hooft coupling which we will define below equation (12.Ih . 


Thus it does not depend on either p or T if p < 1/T, and the size of the instanton is a 
moduli in this region. 
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However, both perturbative QCD and numerical calculation in lattice gauge theory pre¬ 
dict different results. Perturbative QCD predicts that the instanton density for a small 
instanton at T = 0 is 

f StT^ \ HAT C 

n(p,0) oc exp oc p 3 , (1.2) 

where g^(p) is the coupling at scale p [H]. Thus small instantons are suppressed. At 
high temperature T ^ Tc in the deconfinement phase, because of the electric screening, 
large instantons would be suppressed. Indeed the perturbative calculation shows the large 
instanton suppression as 

n(p,r) = n(p,0)exp - log (^1 + (1.3) 

for 71 pT ^ 1 and g^{T) -C 1 [H]. Although the perturbative calculations are valid only in 
limited parameter regimes, such suppressions of small and large instantons would hold for 
any temperature in the deconfinement phase, and then the instanton density would have a 
peak at a finite value. Actually this tendency has been observed in lattice calculations in 
the deconfinement phase [T3] . 

These results clearly disagree with the holographic result (II.ip in the black D4-brane 
geometry. Although holographic QCD cannot reproduce the actual QCD results quantita¬ 
tively in principle [2], qualitative aspects of QCD are expected to be captured. Hence this 
discrepancy is a serious puzzle in holographic QCD. More recently, it has been argued that 
the black D4-brane geometry cannot be identified with the deconfinement phase in four¬ 
dimensional QCD; rather, a geometry called localized solitonic D3-brane will correspond to 
the deconfinement phase in QCD [1]. In this article, we study the instantons in this new 
setup, and see that the instanton density obtained from the localized D3-brane geometry 
satisfies the expectations from QCD. 

We find that the size distribution of the instantons is peaked at a finite value, and 
becomes delta-functional at large-A^. Interestingly, the topology of the gravity dual of the 
stable instantons differs from that of the small ones. Also, we will see that fluctuations 
of the topological charge, which is large in the confinement phase and suppressed in the 
deconfinement phase, would smoothly change at the Gross-Witten-Wadia type (GWW) point 

[MiiniiiH]. 

This paper is organized as follows. We begin in section |2] by reviewing the holographic 
QCD at finite temperature and discussing the geometries corresponding to the confinement 
and deconfinement phases. Then in section [3l we argue instantons in the confinement geom¬ 
etry. These two sections are mostly a review of the previous studies. In section 01 we argue 
instantons in the deconfinement phase. We also argue the 6 dependence and topological 
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susceptibility in section 0 and show that a continuous transition of the susceptibility occur 
at the GWW point in section [ 6 l 


2 Confinement and deconfinement phase in holographic 
QCD 


In this section we review the conhnement and deconhnement phases in four-dimensional 
SU{N) pure Yang-Mills theory in Witten’s holographic QCD model [2]. Let us hrst consider 
a ten-dimensional Euclidean spacetime, whose a: 4 -direction is compactihed on a circle with 
period L 4 , which we call We consider N D4-branes wrapping on this circle. For the 
fermions on the branes, we take the anti-periodic boundary condition along so that 
supersymmetry is broken. 

By taking the large N limit of this system a la Maldacena at low temperature [I9l|20], we 
obtain the dual gravity description of the compactihed hve-dimensional SYM theory on the 
D4-branes [2], which consists, at low temperature, of a solitonic D4-brane solution wrapping 
the Sl^. The explicit metric and dilaton is given by [20] 


ds^ = a' 


u 


3/2 




TT 


dt^ -I- E dxl f4{u)d: 


/,(„) = 1 _ 


2=1 







-I- v?dVtl 


{2tt)^N y ^Ag/dTT 


( 2 . 1 ) 


This solution also has a non-trivial hve form potential which we do not show explicitly. Here 
Ap+i is the’t Hooft coupling on the Dp-brane world-volume, which is given in terms of the 
string coupling gg and Regge parameter a' as Ap+i = . We will also use 

the dimensionless coupling Xym = 2 A 5 /L 4 . 

Since the a; 4 -cycle shrinks to zero aX u = Uq/vo. order to avoid possible conical singularities 
we must choose the asymptotic periodicity L 4 as 


L4 ^/ XYf ^ -1/2 

27r 3 ° 


( 2 . 2 ) 


With this choice, the contractible a: 4 -cycle, together with the radial direction u, forms a 
so-called cigar geometry, which is topologically a disc. Note that this gravity solution is 
reliable in the regime \ym S> 1 where the stringy corrections are suppressed. 

Witten pointed out that four-dimensional pure Yang-Mills theory is obtained in a regime 
Xym -C 1, because the KK modes about and matter helds (fermions and adjoint scalars 
which acquire masses via loop corrections) in the hve-dimensional super Yang-Mills theory 
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are decoupled. Although this QCD regime {Xym ^ 1) and the strong coupling regime 
{Xym 3> 1), where the gravity analysis is reliable, are completely opposite, their properties 
would be qualitatively related as far as no transition occurs between them. (This is analogous 
to the strong coupling expansion of the lattice gauge theory.) Indeed there is a lot of evidence 
which supports this connection, and we expect that supergravity analyses capture qualitative 
aspects of large-QCD. 

So far, we have considered four-dimensional pure Yang-Mills theory at zero temperature. 
In order to study properties at finite temperature, we compactify the Euclidean temporal 
dimension to a circle, and identify its circumference (3 with the inverse temperature, (3 = 
1/T. In four-dimensional theories with fermions, the anti-periodic boundary condition along 
the temporal circle is imposed for fermions. In Witten’s setup, however, fermions in five¬ 
dimensional theories decouple in the four-dimensional limit Xym —t 0 , and hence we do 
not have to impose the anti-periodic boundary condition. Rather, ref. [1] argued that the 
periodic boundary condition is more useful to investigate the QCD deconfinement phase via 
super gravity. 

As we decrease f3 = 1/T in the geometry (12.ip . it reaches O {L4^/ \/Xym) , which is the 
order of the effective string length at u = mq M. Below this value, winding modes of the 
string wrapping on the /3-cycle could be excited. Thus the gravity description given by fl 2 .ip 
would be valid only if 


/ 5 » 


u 

VXym 


(2.3) 


In order to avoid this problem, we perform the T-duality transformation along the f-cycle 
and go to the IIB frame, where the solitonic D4 solution becomes solitonic D3-brane solution 
uniformly smeared on the dual t-cycle0 From now, t' and /3' denote the dual temporal 
coordinate and its period 


13' = = {2'KfT. (2.4) 

^Note that the T-duality along the t-cycle maps the system to the IIB string theory since the periodicity 
of the fermions along this cycle is taken to be periodic. If we took the anti-periodic boundary condition, 
the system is mapped to the OB string theory in which the brane solution has not been studied well. Hence 
Ref. [T] took the periodic boundary condition and focused on the IIB supergravity. However, it may be 
possible to derive similar results from the OB theory too. Another difference in the case with the anti- 
periodic boundary condition is the existence of the black D4-brane solution which is stable for T > I/L 4 
at strong coupling {Xym ^ 3/Li). (Note that the black D4-brane solution is not allowed if we took the 
periodic condition.) Although this solution is thermodynamically favoured at strong coupling, this solution 
is not related to the four-dimensional QCD [UHS]. (Roughly speaking, this solution is an analogue of the 
“doubler” in the lattice gauge theory.) Indeed this solution is not stable in the weak coupling where we 
obtain the QCD [TJ [15], and we should remove this solution by hand if we study the QCD through the 
holographic QCD with the anti-periodic boundary condition. 
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Figure 1: The schematic relations among the free energies of various solutions in the GL 
transition and their topologies. The red lines denote the stable solutions. The blue lines 
denote the meta-stable solutions. The yellow lines denote the unstable solutions. The 
solitonic D4-brane (uniformly smeared solitonic D3) becomes unstable at the GL instability 
point. The non-uniformly smeared solitonic D3-brane solution appears at this GL instability 
point and merges to the localized solitonic D3-brane solution at the merger point. 


In this frame, the mass of the winding strings become heavier as /5 decreases (hence the dual 
radius jS' increases) as opposed to those in the IIA frame, and we can explore the model at 
higher temperature. 

As 13 decreases, the uniformly smeared solitonic D3-brane solution becomes thermody¬ 
namically unstable at a certain temperature, which is called the Gregory-Laflamme (GL) 
instability point [21] p and is numerically given by 


(3GL,inst — 14.4 


A 


YM 


(2.5) 


See also Figured] It is expected that, even before (3 is lowered all the way down to (3GL,inst, 
the smeared solitonic D3-brane solution becomes meta-stable and undergoes a hrst order 
Gregory-Laflamme (GL) transition at an inverse temperature Pgl which is approximately 
given by [1] 


/3gl ~ 



17.1- 


'VM 


(2.6) 


leading to a more stable configuration of D3-branes localized on the dual cycle, whose topol¬ 
ogy is different from the smeared D3-brane solution. See Figure [TJ This hgure also shows 
that the localized solitonic D3-brane solution ceases to exist if 13 is too large; intuitively, if 
f3 were too large, the dual cycle would become smaller than the size of the localized soliton, 

^ GL instabilities have been studying in black strings which are the double Wick rotation [t' 0 : 4 ) of the 

smeared soliton. As far as thermodynamical properties, we can read off the soliton results from the black 
hole ones. 
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which is not possible. The metric of the localized solitonic D3-brane geometry is approxi¬ 
mately given by that of D3-branes on x for a sufficiently large radius f3' of the dual 
cycle [22l [23], which we will see in section 14.11 

It is argued in [1] that this localized solitonic D3-brane geometry can naturally be re¬ 
garded as a counterpart of the deconfinement phase in QCD, and the confinement/deconfinement 
transition can be identified with the GL transition with the transition temperatur^ 


T = 


I^GL 




(2.7) 


In addition to the uniformly smeared D3-brane and localized D3-brane solutions, there 
is another solution: solitonic D3-brane non-uniformly smeared on the f'-cycle. This solution 
describes the D3-branes localized on the f'-circle but there is no gap. (See Figured]) Thus 
this has the same topology to the uniformly smeared solitonic D3-brane geometry while the 
translation symmetry is broken. The metric of this solution is perturbatively derived in |22] . 
This non-uniform solution arises at the GL instability point. Although the behavior of the 
non-uniform solution has not fully understood, it is expected that this solution merges with 
the localized solitonic D3-bran solution as shown in Figured] [25] [33] jSU [35] [36]. This point 
is called “merger point.” 

What is the corresponding phase to the non-uniform solution in QGD? Recall that we 
have taken the T-duality along t-circle and the T-duality maps the locations of the branes 
to the eigenvalues of the Polyakov loop operator exp Aq j. Thus the non-uniform D3- 
brane solution describes a phase characterized by the non-uniform eigenvalue distribution of 
the Polyakov loop operator. Indeed such a phase is well known in large-gauge theories 
although it may be unstable [ST] [38] [39] . In particular, the merger point is an analogue of 
the GWW point in the two-dimensional Lattice gauge theory [161 [HI [H] • In section [6] we 
will discuss the importance of the merger point for understanding how the difference of the 
topological fluctuations at low and high temperatures arises. 


3 Instanton in confinement phase 

Now we consider instantons. First we review the instantons in confinement phase in holo¬ 
graphic QGD. In the bulk theory, the QGD instanton corresponds to the DO-brane winding 

^Indeed there are various evidences which show the resemblance between the GL transition and the 
confinement/deconfinement transition. We can show that the phase transition in the five-dimensional SYM 
theory at strong coupling indeed occurs around (12.61) by applying the analysis in [23] . Several calculations in 
low dimensional gauge theories also agree with this proposal [25] EU ET] El] . Besides |29l EOl EH [32] revealed 
that the confinement/deconfinement transitions exhibit similar properties to the GL transitions. 
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DO-brane 



a) DO-brane in the solitonic D4 geometry b) DO-brane in the black D4 geometry 


Figure 2: Brane configurations of the QCD instanton: (a) DO-brane in the solitonic D4-brane 
geometry, which describes the conhnement phase, and (b) DO-brane in the black D4-brane 
geometry, although it does not correspond to any QCD phases. In the conhnement case, the 
DO-brane can shrink to a point at the tip of the soliton and the DBI action becomes zero. 
In the black D4-brane background, the effective length of the DO-brane does not depend on 
the position, and hence, the potential is constant as in equation fll.ip . 


on the X 4 -circle [9l[ini[II]. Then the brane conhguration of this system is summarized as 

(0) 1 2 3 (4) 5 6 7 8 9 

N D4-branes _____ 

DO-brane (QCD instanton) — 


Here the parentheses denote the compact directions. To investigate the potential for an in¬ 
stanton we evaluate the DBI action of the single DO-brane. Ref. [T0l[TT| showed that the DBI 
action in the solitonic D4-brane geometry fl2.ip . which will correspond to the conhnement 
phase, is 




A 


YM 


Mn 




(3.2) 


where u is the position of the DO-brane along the radial coordinate. Therefore the DO-brane 
is attracted toward the horizon of the D4-soliton [u = Uq), and the classical action disappears 
when it arrives at the horizon]^ This result can intuitively be understood through the cartoon 
of the brane conhguration depicted in Figure [2] (a). Since the DO-brane can shrink to a point 
at the tip of the soliton, the DBI action becomes zero there. 

Now we interpret this result as the corresponding QCD instanton dynamics. The position 
u of the DO-brane would be related to the size p of the instanton mm- In the case of the 

^If we evaluate the backreaction of the DO-brane, we see that the energy of the DO-brane is not exactly 
zero. The solitonic D4-brane solution with non-zero DO-brane charge (D0-D4 geometry) has been calculated 
in [To], and the DO-branes cost the energy ~ VaLlng/Ag, where no = Nq/PVo is the charge density for 
DO-branes and V 3 is the spatial volume of in QCD. Here we have assumed that the density uq is small 
and uniform on the four-dimensional space. See also |40l |41] for the application of the D0-D4 geometry to 
holographic QCD. 
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extremal Dp-brane geometries, this relation can be understood explicitly. Roughly speaking, 
in these geometries, the typical energy scale at the radial position u is given by 
where Xp is the ’tHooft coupling on the Dp-branes [201 021 03] and hence the size of the 
instanton is related to its inverse . Indeed this relation has been conhrmed 

in the AdS 5 /CFT 4 case [H HU HU 06l However this argument cannot be applied to 
the conhnement geometry fl2.R at hnite temperature, since there are two energy scales 
^u/L 4 ^Xym and ^/uqJTmX^u (Through fl2.2p . the latter becomes \/uo/L/^X ym ~ I/T 4 , 
which is the same order to the glueball masses in the holographic QCD [3].) Although we do 
not have explicit relation between p and u in the confinement geometry, we naively assume 

mm 


L4^Xym Q\ 

P~\/ -• 3.3 

V u 

This assumption would be valid at least when Uq where Uq would be irrelevant or when 
M ~ Mo where the two energy scales are coincident. Once we admit this assumption, the 
equation (13.2p indicates a suppression of a small instanton in the conhnement phase, which 
is qualitatively consistent with the perturbative QCD ( 11 . 2 p . At m ~ Mq (p ~ \/LaXymT^ ~ 
L4), the instanton can exist with the zero value of the DBI action, which would imply that 
the huctuation of the topological charge is large [9]. Moreover a larger instanton cannot 
exist. Thus the instanton density has a sharp peak at 


Ppeak ~ T 4 (3-4) 

for large N. The location of the peak does not depend on temperature. This would be 
related to the large-A^ volume independence [50l 1^ . Remarkably the lattice calculation 
in the conhnement phase yields a similar sharp and temperature independent peak in the 
instanton density [I3] . 


4 Instanton in deconfinement phase 

To investigate the thermodynamics of QCD at high temperature through holography, we 
need to take the T-dual along the Euclidean time circle as argued in Section [2l Then the 

® In the AdS 5 /CFT 4 correspondence |48l|49], we read off (TrFAF) in the gauge theory from the boundary 
value of the RR scalar field which is sourced by a D-instanton in the bulk HU HU HSl SI] ■ Hence the wave 
equation of the scalar in the AdSs fixes the instanton dynamics on the boundary. Importantly the wave 
equation for the S-wave can be rescaled so that it is described by a single dimensionless parameter X^k'^/v? 
where k is the longitudinal momentum [43]. It leads us to the scaling behavior p ^ I/fc ~ y/X^/u of the size 
p of the instanton corresponding to the D-instanton located at u. Similar scaling, with the dimensionless 
parameter Xpk^/u^^~P\ would be obtained in other extremal Dp-brane cases too. However the metric of 
the confinement geometry m involves the additional factor / = 1 — (uq/u)^ and the scaled wave equation 
depends on both X^kX /u and X^kX/ uq. Thus we obtain two energy scales y/uJTX^XYM and ^ u^j L/^Xy m- 


















a) Dl-brane in the smeared D3 geometry b) Dl-brane in the localized D3 geometry 


Figure 3: Brane configurations of the QCD instanton in the T-dual picture: (a) the conhne- 
ment phase (Dl-brane in the smeared D3-brane geometry) and (b) the deconhnement phase 
(Dl-brane in the localized D3-brane geometry). In the conhnement case, the Dl-brane can 
shrink as the DO-brane in the solitonic D4-brane geometry. In the deconhnement case, the 
Dl-brane remains a hnite volume and the DBI action is always non-zero. 


brane conhguration (13.Ih is mapped to the IIB frame: 

(O') 1 2 3 (4) 5 6 7 8 9 

N D3-branes _ _ _ _ (^ 4 _X) 

Dl-brane (QCD instanton) — — 

Thus we should consider a Dl-brane instead of a DO-brane to study the dynamics of the 
QCD instanton. 

At T < Tc, the stable geometry in the IIB frame is the uniformly smeared D3-branes, 
which is the T-dual of the solitonic D4-brane geometry (12.Ih . and the instanton is described 
by a Dl-brane on this geometry. (See Fig. |3] (a).) Since the T-duality retains the values of 
the classical action, the results in the conhnement phase discussed in the previous section 
remain the same. 

At T > Tc, the stable geometry is given by the localized D3-branes. We study the 
dynamics of the Dl-brane on this geometry by using the probe approximation. 

4.1 Geometry of localized D3-branes on a circle 

First we explain the details of the localized solitonic D3-brane solution, which corresponds 
to the deconhnement phase. In this geometry, the D3-branes are localized on f'-cycle, where 
t' is the euclidean time direction in the IIB frame. We set the location of their center of mass 
to be t' = 0. Then because of the periodicity t' = t' + [3' {f3' = dvr^T), their mirror images sit 
around t' = n(3' {n = ±1, ±2, • • ■). Since the D3-branes and their mirrors are gravitationally 
interacting, each “horizon” is stretched along t'-direction, and the spherical symmetry is 
broken. Because it is difficult to treat this ehect exactly, we use an approximation which is 
justihed at T/Tgl 3> 1, where the interaction becomes weak. In particular, at the leading 
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order of this expansion, we can treat the horizon spherically symmetric, and we do not 
consider higher order corrections in this paper for simplicity]® 

Althongh snch localized solitonic solutions have not been investigated well, localized black 
brane geometries, which are just the double Wick rotation of the solitonic brane geometries, 
have been studied very actively in the context of the Gregory-Laflamme instability, and we 
can borrow the results. We consider the wick rotation of the black D3-brane localized on 
S'^-circle [23] , and then, the metric takes the following form 


ds^ = 


dxl + fdx 


i=0 


^1/2 




A 


Kd-‘ 


:dv‘^ + KR^dnl 


where 


/ = i 


IP' 


(4.2) 


(4.3) 


and A and K are functions of R and v. Here, X 4 is not the euclidean time direction, but that 
for Sl^ as in the IIA frame. The euclidean time coordinate, t', is included in the (R, n)-plane. 

In order to sinmlify the analysis, we consider two regions: asymptotic region and near 
region [22l [231 [52] h| In the asymptotic region (i.e. u S> or t' 3> uh), effects of the black 
hole can be calculated by solving linearized equations, and then, the metric is given by 


ds^ = a' 


H ^ dx\ + (1 + 2^)dxl I + — -<h) [du^ + + u^dfll) 


. i=l 




2 A 5 / 


+ (d — 


= 


2tiN(5' 


$ = 




+ it' — nl3'y 


Uh 


— \/ 2 A^ 




Note that mirrors contribute to the metric. In the near region (i.e. uR' ~ uh)-, the effect of 
the black hole becomes much larger than that of the mirror images. In this limit, the metric 
is given by 


ds^ = ol 




H = 


2A5 /fd 


2=1 




-2 


1 + 5 


dxl 


To = 


Uh 

V 2 


TT 


+ + (dr^ + RdQl) 




Ar 


2 L 4 


= 


27rAr/3’ 


(4.5) 


® Although this approximation is valid only for high temperature T/T^ ^ 1, the qualitative properties of 
the localized D3-brane would not be changed even around R as indicated in the numerical calculation of 
the localized black holes |55] . 

^ In this paper, we use the approximated form for localized neutral black holes in |52| . It is straightforward 
to obtain the geometry for D3-branes from that for the neutral black holes [22] [23] . 
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where we have dehned the coordinate r such that it approaches to the Newtonian gauge as 
r » Tq. 

In the following, we introduce a midpoint ri, and treat r < ri and r > ri as near and 
asymptotic regions, respectively. 


4.2 The Dl-brane in the asymptotic region 


We hrst consider a Dl-brane located in the asymptotic region where the background metric 
is approximated by fl4.4l) . By regarding the brane configuration fl4.1l) and the symmetry of 
the background geometry, the Dl-brane will be embedded in t', X 4 and u space. We take 
{t', X 4 ) as the world volume coordinates on the Dl-brane, and then the induced metric is 
given by 




+ 24)(ix| + 




(4.6) 


Then the DBI action in the asymptotic region becomes 




fL4 


^Dl = 


(27r)Q!' 

NPL4 

Ak 


dt' / dx^e det gm, 




dt' [1 + 


fWf 


{u:^ (]' :$> uh) (4.7) 


Since t' dependence of $ can be neglected for large u, we approximate that U{t') is constant. 
Then we obtain 


_ Stt^A^ _ XymT 
Xym 2® 


(4.8) 


The second term indicates that the Dl-brane is attracted toward the D3-branes (n = 0). This 
is similar to the conhnement geometry case fl3.2l) . although the potential is now proportional 
to temperature. Around T ~ Tc, this term becomes ~ Nu^/Xymu'^ which is the same order 
to the attractive potential in the conhnement phase fl3.2p at large n, and it becomes stronger 
as temperature increases. 

This result is valid only for u S> uh and the approximation becomes worse as u approaches 
to Uh- At u Uh, 04.8p behaves as 


_ Sn'^N _ N [¥, 
Xym Xym V T 


and for u < uh the above discussion will completely break down. 


(4.9) 
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4.3 The Dl-brane near D3-branes 


Since the Dl-brane is attracted toward the D3-branes, the Dl-brane would be stabilized at 
M = 0 and would stretch between the D3-branes and their mirror image along the compact t' 


circle as depicted in (b) of Fig. [3l In Appendix A we demonstrate that the stable classical 
solution of the DBI action is given by this conhguration indeed. 

Note that the hypersurface at m = 0 of the localized D3-brane geometry has a topology 
of X E? X and the stable Dl-brane wraps on this On the other hand, the Dl-brane 
in the asymptotic region {u S> uh) winds t'- and X 4 -cycles which compose a topology of 
T^. Thus the topology of this Dl-brane differs from that of the stable Dl-brane at m = 0 
which winds This is because the stable Dl-brane reaches the “horizon” of the D3-brane 
where the X 4 -direction shrinks. It indicates that the Dl-brane in the asymptotic region 
cannot continue to the stable Dl-brane at m = 0 smoothly. When the Dl-brane reaches 
the “horizon,” the topology changes. We will later see that the Dl-brane in the asymptotic 
region describes a small instanton, and it means that the small instanton does not smoothly 
continue to the stable instanton. 


Let us compute the value of the DBI action for the stable Dl-brane at m = Olj Near the 
D3-branes or their mirror image, the metric can be approximated by that for the near region 
fl4.5p . However, around the middle between the D3-branes and their mirror, the metric 
cannot be described by that for the near region but should be approximated by that for the 
asymptotic region fl4.4p . 

To evaluate the DBI action for the stable Dl-brane at m = 0, it is convenient to rewrite 
the metric fl4.4l) and fl4.5p in the following combined expression; 


ds^ = a' 


H dxj + fidx\ I -|- idr'^ + r^dQl) 


. ^=1 


(4.10) 


where fr and H are functions of r and one of the angular coordinates of S^, which are 
related to u and t'. In the asymptotic region, they approach to 

2A5//3 


/4 — 1 + 2 $ , 




^ = 5 : 


-I- (F — njd'YY ’ 


(4.11) 


with = V? E In the near region, the geometry has spherical symmetry on at the 

®If we could calculate the value of the DBI action for the appropriate configuration of the Dl-brane 
corresponding to the QCD instanton with size p, we would obtain the potential for p as we did for the 
solitonic D4-brane background in section [31 Since the localized D3-brane does not have the isometry along 
t'-cycle, it is difficult to specify the configuration for a specific size of the instanton. For this reason, we 
calculate the action only for the stable classical solution. We will discuss a related issue in section IT4l 
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leading order and they become 

/ 1 _ ^ ^ 
U = 


' 0 
7^ 


1 + ^ 


/. = l + i, H = 


2 A 5//3 


^4\ -2 


(4.12) 


as is shown in fl4.4p and fl4.5p . 

In this metric, u = 0 corresponds to a hxed direction in and r can be chosen to be 
identihed to t' when u = 0. Then, the induced metric on the Dl-brane at m = 0 is expressed 
as 

= a' f^dxl + fr{t')dt''^~\ , (4-13) 

and the DBI action is given by 


S'di — 2 


N(3 


r-/372 


-.L 4 


dt’ I dx^^/JJ'r 


(4.14) 


'I’D 


In order to calculate this action, we introduce a mid-point ri and divide the t'-integration 
into two parts, that for the asymptotic region and that for the near region; 


•S'di — iS'near + S-, 


asymp • 


(4,15) 


We will soon see that the hnal result is independent of ri. 
The integration for the near region can be calculated as 


7 — 9 


= 2 


A 5 

Ar 


dt' 


NpL, 1-^ 


I'D 


1 10 . 
i + i/4 


V2 + ^ + 0( 


ro 


ro 


ri 


(4.16) 


Here, we are assuming tq <C 13', and hence, we can take ri ^ tq. In the asymptotic region, 
we can neglect <h since it gives contributions at 0 (^ 0 ), and hence the DBI action can be 
calculated as 




2N/3L, 


NI3U 


asymp 




dt' + O(rii) = (/?' - 2n + 0(ri)) 


'ri 


A,r 


By summing these two results, we obtain 

7V/3L4 


/S'di — 


Ak 


{I}' -2'''^r„ + 0(rl)) . 


(4.17) 


(4,18) 


which does not depend on ri. By using f3' = (27r)^T, we hnally obtain 

Svr^iV 


/S'di — 


A 


1 - 


YM 


^YM 

dvr V T/f 


Stt^N 

^YM 


1 - 


3^/2 /t 


2iV2vrV T 


(4.19) 
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Thus the DBI action is hnite and the topological fluctuation is exponentially suppressed at 


a 


large-iV. 

Note that the value of the action fl4.19p for the stable Dl-brane at u = 0 is the same 
order to (14. 9 h for u ~ uh which is extrapolated from the DBI action (14.Sh for large u. It 
would indicate that the potential (14. 8 p at large u continues to the value (I4.19p . Recall that 
the topologies of the Dl-brane at large u and u = 0 are different, and the topology change 
occurs when the Dl-brane reaches the “horizon” of the soliton. Since the value of the classical 
action is related to the area of the Dl-brane, it would be continuous through the topology 
change. However, its (higher) derivative with respect to some deformation parameters of the 
Dl-brane may not be continuous. 


4.4 Size of instantons in the deconfinement phase 

We have calculated the DBI action of the Dl-branes. Now we argue the relation between 
the size p of the QCD instantons and the radial location u of the Dl-branes, as we have done 
for the confinement phase in section [3l The relation is more complicated than that for the 
confinement geometry (12.ip . since we have taken the T-dual on the temporal circle and the 
energy in the IIA frame would appear in an unusual manner. Furthermore, the metric can 
analytically be expressed only by a couple of the approximated forms for two patches. 

Fortunately the asymptotic metric fl4.4p has an approximate isometry along the temporal 
circle if u is sufficiently large, and we can consider the IIA frame by taking the T-dual again. 
There the typical enegy scale is L 4 XYM for the Dl-brane which is located at u. This is 
the same scale to that of the confinement geometry, since the localised D3-brane geometry 
fl4.4p asymptotically approaches to the smeared D3-brane geometry [T], which is the T-dual 
of the confinement geometry fl 2 .ll) . 

For the near D3-brane metric fl4.5l) . the isometry along the temporal circle is broken. 
Hence the T-dual picture in the HA side is not clear and it is hard to see the relation. If 
we consider only the near region, there typical energy scales would be naively r/ 
and th/\/XymL 4 T ~ I/T 4 for the Dl-brane located at r. However we need to consider 
the connection to the asymptotic region, where the aforementioned different scalings arise, 
to estimate the energy of the gauge theory on the boundary. (See footnote E] and [45]. 1 
Furthermore, the Dl-brane for the stable configuration at n = 0 is stretched along r direction, 
and it makes the situation more complicated. 

However, for small instantons the contribution of the near region would be irrelevant and 
the asymptotic metric (14.4p would dominate for obtaining the relation to the instanton size. 

® If we use the dilute gas approximation, we obtain xt oc in the localized D3-brane geometry where 

Xt is the topological susceptibility. 
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Then it would be possible that the p — u relation fl3.3l) in the confinement geometry holds 
approximately even in the deconfinement phase. Under this assumption, the DBI action 
fl4.8p is rewritten as 


<S'ni — 


Stt^AT 


A 


YM 


'^YM 


+ ••• 


(4.20) 


The value of this action is larger for smaller p and it suggests that the small size instanton 
would be suppressed. 

For the stable Dl-brane at m = 0, we cannot estimate the instanton size because of the 
difficulties mentioned above. However the distance between the Dl-brane at m = 0 and the 
boundary [u = oo) is finite, which implies that the instanton size for the Dl-brane at m = 0 
would be finite. We presume that this corresponds to the largest instanton in QCD and a 
larger instanton is not allowed. Thus the instanton density would have a sharp peak at this 
value of p at large-iV. 

Although we cannot calculate this largest size, we can estimate the lower bound for this 
size by using the relation (13.31) in the asymptotic region and substituting u = uh. 


p ~ 


\YMJ^4 


Uh 



(4.21) 


Note that this lower bound of the peak size (I4.2ip increases as T decrease and reaches p ~ L 4 
around T ^ T-, which is the same order to the peak size of the instanton in the confinement 
phase (131 


5 ^-vacuum and topological susceptibility 


We have studied the instantons in the deconfinement phase. The results show that the 
instanton density has the sharp peak at a finite instanton size but the energy at this size is still 
finite. This implies that the topological fluctuation would be suppressed in the deconfinement 
phase. On the other hand, the zero energy of the instantons in the confinement phase implies 
the large fluctuation of the topological charge. To confirm this picture, we investigate the 
instanton effects in 0 -vacuum and estimate the topological susceptibility Xti which is defined 
by the second derivative of the free energy with respect to 9 parameter: 


_ (fF 

“ IF ' 


(5.1) 


Recall that the value of the DBI action at u ^ uh (14.91) obtained from the potential (14.81) for large u 
is the same order to the DBI action at u = 0 (j4.19l) . Thus, the size (I4.2ip obtained from the relation (13.31) 
at It ~ Uh provides the order of the largest size of the instanton (the Dl-brane at u = 0), if it has a similar 
property. 
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We show that the topological susceptibility is indeed suppressed in the localized D3-brane 
geometry consistently with the hnite value of the DBI actioiJ^. 

We hrst recall the topological charge at low temperature mm- The conhnement phase 
corresponds to the solitonic D4-brane geometry, and the instanton is described by the (eu¬ 
clidean) DO-brane wrapped on the X 4 -direction. Thus the topological charge corresponds to 
the “RR-charge” and can be estimated from the conhguration of the RR 1-form Ci. The 
parameter 6 , which is the chemical potential of the instanton, corresponds to the boundary 
condition of Ci; 

e = [ Cl, (5.2) 

where integration is over which is a; 4 -direction at the boundary u —)■ oo. In the case 
of the solitonic D4-brane geometry, is the boundary of a disk D. By using the Stokes 
theorem, it can be written in terms of the held strength F 2 = dCi; 


e = [ F 2 . (5.3) 

J D 

This implies that the held strength is proportional to 9. Then, the classical action for the 
bulk RR-held is estimated as; 

^ ~ j F^(xe\ (5.4) 

Therefore, the free energy has hnite quadratic term of 6 , indicating the topological suscep¬ 
tibility Xt is hnite. Thus the huctuation of the topological charge is large [9]. 

Now, we turn to the localized D3-brane geometry. In this case, the instanton corresponds 
to the Dl-brane which is wrapped on the torus of the ( 0 : 4 , f')-plane. The parameter 9 is 
related to the boundary condition of the RR 2-form C 2 ] 


9 = [ C 2 (5.5) 

where the integration is performed at the boundary u —)■ 00 . In order to see the contributions 
of this boundary condition to the free energy, we consider the Stokes theorem; 



(5.6) 


where T 3 is the held strength associated to C 2 and M is the three-dimensional space of 
{u,X 4 ,t'). The circle of X 4 shrinks to a point at r = tq, which can approximately be 
expressed in the (n, f')-coordinates as 


+ 1 '^ 


^^For lattice studies, see e.g. [Ml 1551 [5S] . 


(5.7) 
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Thus the X 4 -direction can shrink only at \t'\ ^ tq. The space continues to u = 0 in the other 
region, and connected to the opposite side of Then, the boundary of M consists of two 
tori, and at u = oo with opposite angles, and the Stokes theorem provides us with 



(5.8) 


where the minus sign in the last term comes from the difference of the orientation for two 
tori. The 2-form held C 2 takes the same value on and since they are the opposite 
points on and we assume that the solution has the spherical symmetry. This implies 
cancellation in the r. h. s., and hence, the held strength is not constrained by 9. Thus the 
topological susceptibility Xt vanishes and the huctuation is suppressed as we expected!^ 


6 Continuous transition of topological fluctuation at 
GWW point 

As we have shown in section 01 the DBI action of the Dl-brane in the localized solitonic 
D3-brane geometry remains hnite and the topological charge huctuation is exponentially 
suppressed even at the critical temperature (12.61) . This is not surprising since the conhne- 
ment/deconhnement transition (GL transition) is of hrst order. 

As in Fig. [H the localized D3-brane branch is connected to the conhnement geometry 
(uniformly smeared D3-brane) through the non-uniformly smeared D3-brane geometry. By 
tracking this, we can see how the physics in the deconhnement phase changes to that in the 
conhnement phase. Then an important question is where and how the suppression of the 
topological huctuation in the localized D3-brane geometry changes to the large huctuation 
in the conhnement geometry. We propose that it will occur at the merger point where the 
localized solitonic D3-brane geometry merges to the non-uniformly smeared solitonic D3- 
brane geometry and the topology changes. Since there is no “gap” along t'-circle in the 
non-uniform D3-brane geometry similar to the uniform D3-brane geometry, the DBI action 

The black D4-brane geometry also shows Xt = 0 nn, even though it does not provide the correct 
instanton density. In this case, (u,X 4 ) plane is terminated at the horizon, and hence it has the topology of 
the cylinder. The parameter 0 is given in terms of the RR 1-form by 

0= f F 2 - f Cl , (5.9) 

Jc Js\j 

where is the X 4 circle at the horizon. Now F 2 can be zero by adjusting the second term according to 0, 
and then, the topological susceptibility Xt is zero. We can explain Xt = 0 in a similar fashion even in the 
case of the localized D3-branes. If we restrict {u,X 4 ,t') plane to m > 0 by using the spherical symmetry on 
5'^, the integration of C 2 at m = 0 appears instead of the last term in (15.81) . It can be adjusted such that it 
cancels the integration of C 2 at the boundary. 
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of the Dl-brane in this geometry can be zero. On the other hand in the localized D3-brane 
geometry, due to the existence of the gap, the DBI action is hnite. As this gap is becoming 
smaller, the DBI action will be smaller and would reach zero at the merger point. Therefore 
the continuous transition would occur at the merger point. 

This is also consistent with the analysis in the 0-vacuum. If the geometry has only 
single boundary, the held strength of the RR-held is constrained by 9 and the topological 
susceptibility becomes hnite. In the localized D3-brane geometry, a; 4 -direction does not 
shrink in a specihc region and then, the 3-form hux reach to the opposite side of S‘^. This 
ehectively plays the role of diherent two boundaries. However, the “gap” would close at the 
merger point, and then, the hux cannot pass to the opposite side. This implies that the 
topological huctuation is not suppressed at the merger point 

Recall that the merger point will correspond to the Gross-Witten-Wadia type transition 
point in the gauge theory [I61[I71[I8], where the topology of the eigenvalue distribution of the 
Polyakov loop operator changes. It indicates that the topology of the eigenvalue distribution 
of the Polyakov loop is crucial for the topological huctuation in QCD. It sounds reasonable 
since both the Polyakov loop and instantons are related to the conhgurations of the gauge 
helds. 

On the other hand, the translation symmetry along t'-circle, which is broken in both 
the localized D3-brane and the non-uniformly smeared D3-brane geometry, is not critical 
for the topological huctuation. This translation symmetry correspond to the Zn symmetry, 
which characterizes the conhnement [1], and we predict that this symmetry is not directly 
connected to the large huctuation of the topological charge. 
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Appendix A Stable configuration of Dl-brane in the 

localized D3-brane geometry 


In section 14.31 we calculated the DBI action for the stable conhguration of the Dl-brane, in 
which the Dl-brane is stretched between the D3-branes and their mirror images at m = 0. 
Here, we argue that this conhguration is the only stable conhguration of the Dl-brane which 
wraps on t'- and a: 4 -directions. 

We consider the Dl-brane which is embedded in the three-dimensional space of {u, t', X 4 ). 
To investigate the Dl-brane in this space, we express the background metric fl4.10p as 


ds^ = a' 


H I dxl + f^dxl I -|- {dr'^ -|- r'^dO'^ -|- sin^ Odfl^) 

. *=i 


(A.l) 


and choose the coordinate 9 so that the three-dimensional space of {u, t', X 4 ) is parameterized 
by (r, 0, X 4 ), and 6 = 0 (and 6 = ir) corresponds to m = 0. Then the Dl-brane lies on (r, 6 , X 4 ) 
space and the induced metric on it is given by 


ds"^ = OL [H ^^"^f^dxl + (1 + , (A.2) 


where d is a function of r and 6 *' = dO/dr. The DBI action can be expressed as 


Rni — 2 


A 5 


r/372 


-L4 


dr / dx 4 \/f 4 fr{l + r‘^e'‘^) 


(A.3) 


Since 6 dependence of and fr are negligible if tq fd' as we argued in section 14.31 we can 
solve the equation of motion for 9 as 


9' 


Cl 

r\/r‘^Ufr - cj ’ 


(A.4) 


where ci is a constant. 

If Cl = 0, 6 * becomes a constant, and the solution describes the Dl-brane which is or¬ 
thogonal to the “horizon” at r = tq and extends straightly to outside with a hxed angle 6 * 0 . 


19 












Figure 4: Schematic plots of the Dl-brane in {u,t') space. The red lines describe the Dl- 
brane which obeys the equation of motion flA.4j) . 


(See Fig. HI) If we choose 9 = 0 (and tt) so that the Dl-brane lies along u = 0, we obtain 
the stable solution which we investigated in section 14.31 

If Cl 7 ^ 0, the brane is curved in the (r, 6 ')-pIane. For small ci (<^ vq), the solution behaves 
around r = tq as 

^ = ^0 + {Vr - Tc-l -), ^1^0 = ro +-^^ci + 0{cl) (A.5) 

This solution describes the Dl-brane which does not reach to r = tq, but turns at r = Cc 
and goes back to the outside of the near region. (See Fig. |H) Since flA.4p indicates 0' —)■ 0 
as r —)■ oo, the Dl-brane asymptotically extends to angles 6^0 ± 6 *i as r —)■ oo where the 
asymptotic value 6i is fixed by ci. The solution flA.SD indicates that 6i will decrease as Tc 
decreases and achieves 6 *i = 0 at Cc = tq (ci = 0). Indeed we can conhrm that 9i approaches 
to the maximum 7 r /2 as oo by solving (IA.4D explicitly. 

So far we have not considered the periodicity of f'-cycle, and now we impose it to the 
solutions. We demand that the solutions are smoothly connected at t' = ±/9Y2 which are 
the middle points between the D3-brane and its mirrors. This leads the following boundary 
conditions; 


D(F = ±/3'/2) = Uo , 


dU{t’) 


dt' 


t'=±/3V2 


0 , 


(A. 6 ) 


where we have taken the coordinates {u^t') and U{t') is the prohle of the Dl-brane in these 
coordinates. Then we immediately notice that the possible solutions are Ci = 0 with 9 = 0 
and Cl 7 ^ 0 with 9q = 7 r /2 and 9i = ti/2 {tc = oo) only. These are the constant u solutions 
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at M = 0 and u = oo respectively. Thus m = 0 is the only stable solution of the Dl-brane in 
the localized D3-brane geometry. 

However, the higher order corrections of tq might be relevant in the intermediate region 
between the asymptotic region and near region. Although the Dl-brane is approximated 
by straight conhguration in the asymptotic region in the above analysis, the higher order 
correction may bend the Dl-brane and it might allow other solutions. In order to be a 
solution which satisfies the boundary condition flA. 6 j) with Uq ^ 0, it must go toward the 
D3-brane from t' = (3'/2. Let us see whether it happens. The DBI action of the Dl-brane 
in the asymptotic region is given by fl4.7l) . By assuming U' 1, the equation of motion 
becomes 

(A.7) 


d^U{t') 


dt'"^ 


3 ^ 
4 du 


u=U{t') 


We solve this equation around t' = j3 '/2 with the boundary condition flA. 6 j) and obtain 


U{t') ~ U, 


3 ^ 

^ S du 


u=Uo 

t'=0l2 


+ • • • 


(A. 8 ) 


Since > 0 for m 7 ^ 0 from fl4.4l) . the Dl-brane goes away from the D3-branes for Uq 7 ^ 0. 
Therefore the higher order corrections of tq does not change the result. Hence the Dl-brane 
located at m = 0 is the only stable configuration even if we take into account the corrections 
of ro. 

It would be worth comparing with the case of D3-D7 system [53] in which the D7-brane 
has non-trivial stable configurations. In this case, the D7-brane is not straight outside the 
horizon and approaches to U = const. This is because the D7-brane wraps on the and 
hence tends to stay in the region with small radius due to the tension of these directions. 
On the other hand, the Dl-brane does not wrap no cycle other than X 4 and t' and hence 
extends straightly in the near region. 
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